#include <stdio.h>
#include <math.h>

#define SIZE 200
typedef int index;
int main(void) {
	double e[SIZE] = {0.}, h[SIZE] = {0. }, imp0 = 377.0; //Boundary Conditions
	index nt, mx;
       	int maxTime = 250;
	/* Step by step, time step , then x step*/
	for (nt = 0; mx < SIZE - 1; nt++) {
		for (mx = 0; mx < SIZE - 1; mx++)
			h[mx] = h[mx] + (e[mx + 1] - e[mx]) / imp0; 
		for (mx = 1; mx < SIZE - 1; mx++)
			e[mx] = e[mx] + (h[mx] - h[mx - 1]) * imp0;
		e[0] = exp(-(nt - 30.) * (nt - 30.) / 100.);

		printf("%g", e[50]);
	}



	return 0;
}
